Exact mean exit time for surface-mediated diffusion 
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We present an exact expression for the mean exit time through the cap of a confining sphere for 
particles alternating phases of surface and of bulk diffusion. The present approach is based on an 
integral equation which can be solved analytically. In contrast to the statement of Berezhkovskii and 
Barzykin [J. Chem. Phys. 136, 54115 (2012)], we show that the mean exit time can be optimized 
with respect to the desorption rate, under analytically determined criteria. 



I. INTRODUCTION 

Surface-mediated processes, in which a particle ran- 
domly alternates between surface and bulk diffusions, are 
relevant for various chemical and biochemical processes 
such as reactions in porous media or trafficking in living 
cells pn-TTo] . Recently, it has been shown theoretically 
that intermittent dynamics between two diffusive phases 
can lead to faster kinetics. The kinetics was character- 
ized by the mean first passage time (MFPT) of diffusing 
particles to a fixed reactant [TT] . 

As a representative example of confined interfacial ki- 
netics, the case of a particle alternating bulk and surface 
diffusion over a sphere which contains a reactive cap, was 
considered in [T2TU4| . The reactive cap can also be in- 
terpreted as a hole, in which case the MFPT is referred 
to as the mean exit time [T51IT8] . While the desorption 
rate from the surface is independent of any geometrical 
parameter, the switching dynamics from the bulk to the 
surface is determined by the statistics of returns to the 
sphere. These statistics strongly depend on the boundary 
behavior of the process. After each desorption event, it 
was assumed in [TH [T3] that the particle was ejected into 
the bulk at a distance a > that was required to avoid 
an immediate re-adsorption on the perfectly adsorbing 
sphere. A similar description was used in discrete square 
lattice versions of this model, in which the lattice spac- 
ing played the role of the cut-off distance a [HI [20] . In 
these studies, the MFPT was found, under certain con- 
ditions, to be an optimizable function of the desorption 
rate. The favorable effect of desorptions was attributed 
in [12] to the fact that bulk excursions reduce the time 
wasted due to the recurrence of surface Brownian mo- 
tion, by bringing particles through the bulk to unvisited 
regions of the sphere. Previous mean-field treatments 
which ignored spatial correlations, missed this possible 
optimum |21j . In the case of a uniformly semi-reflecting 
sphere, including the target, the MFPT was also found 
to be an optimizable function of the desorption rate even 
for a distance of ejection a set to zero [22] . 

Recently, a coarse-grained approach to the surface- 
mediated search for a perfectly adsorbing target in an 
otherwise semi-reflecting sphere was considered in |23| . 



This model is relevant in numerous real situations in 
which the particle reactivity with the target is not re- 
lated to its affinity with the rest of the surface, as it is 
obviously the case in particular for exit problems. Rely- 
ing on an elegant first order kinetics scheme, concise ap- 
proximate expressions were proposed in |23] for the spa- 
tially averaged MFPTs (called here global MFPTs and 
denoted GMFPTs) for a uniform distribution of starting 
points over either the sphere surface or the bulk. These 
averages will be called hereafter surface and bulk GMF- 
PTs. This approach, which treats the bulk and surface 
diffusive phases as two effective states coupled by first 
order kinetic equations, led to a monotonic GMFPT as a 
function of the desorption rate. This striking difference 
with Ref. [12] was attributed in [53] to the non-locality of 
the desorption process, in which the instantaneous ejec- 
tion at a non zero distance a implied a violation of the 
detailed balance condition. 

In this article we clarify this puzzling situation and ad- 
dress the question of the optimality of the GMFPTs as a 
function of the desorption rate for the mixed boundary 
condition of |23j . More precisely, (i) we provide an ex- 
act solution for the MFPT; (ii) we prove that the surface 
GMFPT can still be optimized with respect to the des- 
orption rate, under analytically determined criteria; (iii) 
we compare our results with the coarse-grained approach 
of [23 , which is shown to be accurate only in a limited 
region of the parameter space. 



II. THE MODEL 

We consider particles diffusing in a three-dimensional 
spherical cavity of radius R (see Fig. [lj switching be- 
tween phases of surface diffusion with diffusion coefficient 
D\ and phases of bulk diffusion with diffusion coefficient 
D 2 - The time spent on the surface is assumed to follow 
an exponential law with desorption rate A. 

The target is a cap defined as the portion of the sphere 
9 £ [0, e], where 9 is the elevation angle in spherical co- 
ordinates. The target is considered as perfectly reactive, 
i.e. particles react as soon as they encounter the tar- 
get for the first time. We will focus here on the mean 
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limits. Indeed a point-like target (e = 0) is detectable 
neither by bulk (3D) excursions nor by surface (2D) dif- 
fusion. 



III. EXACT SOLUTION 



A. Basic equations 



FIG. 1: (Color online). Model of surface-mediated diffusion 
in a sphere of radius R containing a perfectly adsorbing cap (a 
target). For finite values of k, particles (represented as small 
spheres) diffusing in the bulk randomly bind to or bounce 
from the remaining (non-adsorbing) part of the sphere surface 
(outside the target cap). The target is reached either from the 
bulk (shown by the green particle), or from the surface (shown 
by the blue particle). 



first passage time (MFPT) on the target. The surface 
GMFPT (ii) is defined here as the average of the MFPT 
over a uniform distribution of the starting points on the 
sphere surface, including the target. In turn, in Ref. |23j . 
the MFPT was averaged over the starting points out- 
side the target; we denote this average as The 
two averages are related by (ti) = A{ti)[-q, where A is a 
geometrical factor which is written as 



A = 



1 



1 + tan 2 (e/2)' 



(1) 



The non-reactive part of the sphere is considered as 
semi-reflecting: a particle reaching by bulk diffusion the 
non-reactive part of the sphere surface gets randomly ei- 
ther adsorbed or reflected back in the bulk. The proba- 
bility for binding to the surface is an increasing function 
of the adsorption coefficient k, which will be precisely 
defined through the radiative boundary condition for the 
MFPT in Eq. fl5| . In particular, k = oo (resp. k = 0) 
corresponds to a perfectly adsorbing (resp. reflecting) 
boundary. Notice that the model that we consider here 
is exactly the same as in [2"5] . 

The limits A = and A = oo correspond to simple 
situations. In the case A = 0, particles are trapped on 
the surface until they reach the target, and the exact 
expression for the surface GMFPT is [55]: 



t s — (ti)x=o — 
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The limit A = oo is equivalent to a purely reflecting 
boundary. The asymptotic behavior for the narrow es- 
cape limit e <C 1 is given in [24] : 
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Notice that the bulk and surface GMFPTs diverge as e 
tends to zero, as it can be seen in particular in the above 



Using a standard formalism of backward equations 
[27] . we derive the diffusive equations ( 12|13 ) and the 
appropriate boundary conditions ( 14|15 ) satisfied by the 
MFPT, which we proceed to solve in the next section. 

We define p((r, 6, <j>), s\(x', (j)'),t r ) (resp. 
p((0,4>),s\x' ,t')) as the probability for a particle 
located at time t' at x' - with either x' = r', 9' into the 
bulk or x' = 9' on the sphere - to reach at time s > t' 
the point (r, 9, <p) in the bulk (resp. the point (9, <j)) on 
the surface). In what follows, we omit the azimuthal 
coordinates by referring to probabilities averaged over 
the initial and final azimuthal angles: 



p((r, 9),s\x',t' 



#/ d4>'p((r,9,<t>),8\(a>,<l>'),l!) 



o 



o 



For convenience, in this section we will use the shorthand 
notations 

p((r,6),s\x',t') = p(r,d), p(9,s\x' ,t') = P (9). (4) 

Conversely, we define the conditional probabilities: 
p(x,i\(r,9), s) = _p(r,6) and p(x,i\6,s) = p(0), with 
t > s and x — (f , 9) in the bulk or x = 9 on the sphere. 



In the following, the Laplace operator is A( r ,8) 
Ag/r 2 where A r and Ag are 



= A r 



A, 



d 2 



2 d 



1 



dr 2 r dr ' 



dg (sin 9 dg 



(5) 



1. Diffusion equations 

For the process under study the conditional probability 
p(r,9) satisfies the forward diffusion equations, for all 
s>t', [33], 



dp(r,9) 
ds 

d P (9) 
ds 



= D 2 A M) p(r,0) 



R 2 



A e p{9) - A p(0) + kD 2 p(R,i 



(6) 
(7) 



Each term in the right-hand side of Eq. ^ has a straight- 
forward physical interpretation ; these are (from left to 
right) (i) the diffusion within the surface state; (ii) des- 
orption events with a constant rate A; (iii) adsorption 
events with a success rate quantified by k f[22j). 
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Equivalently to Eqs. Q and d7|), the conditional prob- 
abilities satisfy the backward equations, for all t > s' , [27] 

9 J^l = -D 2 A {r , e)P (r,6). (8) 

^ = -^A e p(0) + a{p(0) -p(i2,d)} (9) 

These backward diffusion equations are commonly used 
to determine first passage time observables [28] . 

We define t\{6) as the MFPT for particles started on 
the sphere at the angle 9, and t 2 (r,9) stands for the 
MFPT for particles started at the bulk point (r, 9) (the 
second angular coordinate <f> is irrelevant due to the sym- 
metry and thus ignored). The MFPTs are expressed in 
terms of the conditional probabilities through the rela- 
tions [27] 



h{9) 



t 2 (r,9) 



3G 



dti^J 2nR 2 sin9d9 p(9,t\9,0) 
2nr 2 dr sin 9d9 p(r, 9, t\9, 0) \ , (10) 
dt | / 27ri? 2 sin 6>d0 p(0, t\r, 6, 0) 



ii wo 
2nr 2 dr sin 0d0 p(r, 0, t\r, 8, 0) 
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(11) 



where S = (0,R) x (0,2tt). 

Substituting the latter relations ( 10 1 and (111 into Eqs. 
f[9 ), one can show that the MFPTs satisfy the set of 



equations 
D 



R 2 



j A fl ti(0) + A (t 2 (R, 9) - h{9)) = -1 (e < < tt), 



(12) 



£>2 A 



ta(r,0) = -l ((r,9)eS). 

(13) 

In the next section we specify the appropriate boundary 
conditions for the MFPT. 



2. Boundary conditions 

We justify that the B.C. for the process defined in Sec. 
H] are: 

(i) the Dirichlct boundary condition 



h(6)=0 (O<0<e), 



(14) 



which expresses that the search process is stopped on the 
target. 

(ii) the mixed boundary condition 

f (0 < 8 < e), 

Wo)-\^ <«<»<-0. 

I, k Or \r=(R,e) 



(15) 



in which the first relation expresses perfect adsorption 
on the target (9 £ [0, e\) while the second relation (radia- 
tive B.C.) implements an imperfect adsorption process 
on the rest of the surface (8 <E [e, tt]). The mixed bound- 



ary condition (151 is the major difference of the present 



model from the previously studied case [22] , in which the 
radiative B.C. was imposed over the whole boundary, in- 
cluding the target. Although this modification may seem 
minor, the difference on the target reactivity has a dras- 
tic effect on the surface GMFPT for short adsorption 
times (A » Di/R 2 ) or low adsorption rates (kR <C 1), 
as shown in Fig. [6] below 



In order to justify the form of the radiative B.C. (151, 



we determine the backward B.C. on the probability dis- 
tribution from a well-known forward B.C. For the pro- 
cess defined in Sec. [TT] the forward B.C. equation on the 
probability distribution is [23] 



D 



dp(r,6) 



kD 2 p(R,8) + X P {8). 



(16) 



dr \r=R 

Since the stochastic process under study is Markovian, 
one obtains the following Chapman-Kolmogorov equa- 
tion on the conditional probabilities, for t > s > t' , 

PIT 1>R 

p(x,i\x',t') = / / 2Trr 2 drsin9d9 p{r,9)p{r,9) 



o Jo 



2R 2 KS\n9d8 p(8)p(8). 



(17) 



Taking the derivative with respect to the intermediate 
time s of the above relation leads to the identity 



7T r R 



2,r 2 sin 9d9dr { ^£^p(r, 0) + p(r, 0) ^ ° ] 



V ds 



ds 



2nR 2 sm8d8( d ^lp(9) + d 4 9 lp(e)) = ... 



V ds 



ds 



The next step is to substitute diffusion Eqs. ([6] -[8]) into 
the last relation. The two terms with the angular Laplace 
operators cancel each other due to its hermiticity: 



sin f 



'(A, p(8)p(8)-A e p(0)p(0)) =0. 



The divergence theorem applied on the bulk Laplacian 
A( r ,9) and the backward equations (8j9l yield the follow- 
ing relation over the sphere surface: 

0= r 2R 2 n sin 9d9 ( D 2 9p ^'^ p(r,0) 
Jo V dr l« 

_ D2 dpp0) p(r,9)+\lp(9)-p(R,0))p(d) 
or \r I J 



+ \-\p(8) + kD 2 p(R,8) \p{8) , 



which is satisfied only if: 

D 2 d -^± p(R,8 
or \r 

-Xp(R,0)p(9) 



D- 



dp(r, 8) 



dr \r 
kD 2 p(R, 8)p{9) 



p{R,9) 



(18) 
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Insertion of the forward B.C. ( 16 ) into Eq. ( 18 ) gives the 



B.C. on the backward probability distribution 

dp(x,i\(r,0),s) 



Or 



\r=R 



k{p(x,i\9,s)-p(x,i\(r,0),s)} 



\R- 

(19) 

We then integrate Eq. (|19| over the space and time vari- 



ables x and t according to Eqs. (10) and (11) to obtain 
the B.C. on the MFPT: 

fit 

^ =k{ti(e)-t 2 (R ) 6)} (e<e<7r), 

or |r=(i?,e) 



which identifies with Eq. (15). Agreement with the B.C. 
used in [22] is discussed hiXppendix |A} 



B. Integral equation 



From the set of Eqs. (12) - (15) we now derive an 



integral equation on t\ only. 

We first recall that the eigenf unctions of the angular 
Laplace operator Ag of Eq. |5| are expressed in terms of 
the Legendre polynomials P n : 

- A e P„(cos0) =p„P„(cos0) (n>0), (20) 



where p n = n(n + 1). We set V n (8) = y/2n + lP n (cos9) 
to get the orthonormality: (V^IV^) = S nm , where the 
inner product is defined as 



(f,9)->(f\9h 



f(9)g(9)sm9d9, (21) 



and (f\g) is the scalar product for e = 0. We also define 
K^ n ={V n \V m ) t (m 7 n>0), (22) 
with the explicit expressions listed in Table [I] 

The starting point for solution of the set of equations 



( 12 and ( 13 ) is a Fourier decomposition of t 2 {r, 9), 

+X>»(i?) m*)> ( 23 ) 



h(r,9) = a a 



2dD 9 



where d = 3 for the three-dimensional spherical cav- 
ity considered here (in Appendix [Ej we show how this 
approach can be directly translated for two-dimensional 
problems) . 



Due to the B.C. (15), the projection of t 2 {R,9) onto 
the orthonormal basis {V n (0)} n >o, is 

' t 2 (R, 6)V n {9) sin 6d9 = I t 1 (8)V n {9) sin 9d9 

Je 



k dr \r 



V n (9)sin 9d9. 



(24) 



V n {6) 



F n (u),n> 1 



9e(0) 



g e \l) t = {g e ) e 



Cn (n > 1) 



(<0 



K$ (n>l) 



K { :1 (re > 1) 



Jfil (m / n ) 
n, m > 1 



iki (n > 1) 



Inm (m ^ n) 
m. n > 1 



V2n + 1 P n (cos9) 



E 2(« - l)Pl{u) + [P k (u) - ft_!( U ) 
fc=l 

-(m - 1)P„ 2 M + (it - l)P 2 (n) + it 



1 1— cos(e) 



/2n + l 



In ( - i+f^ 

1 1 — cos e / 2 



p„-l(") 

n+1 



l + cos(e) 

_2_ 



f„ + l(u)-f n -l(») 
2V2ra+l 



F n (u) + 1 
2 



(4"+l)-n7~+l)) « m " n)uP m (u)P n (u) + 

+nP n -i_(u)P m {u) - mP m -\{u)P n {u)) 



2n+l 
2 



»P„(»)-P„-l(.) F„(u)+1 
n+1 ' 2n+l 



2(n + S^(^ + l)T " ™)uP m {u)P n ^ 

+(m+l)P m (u)P n -i(u) 
-(n + l)P»(u)P m -i(u)) 



TABLE I: Summary of the quantities involved in the compu- 
tation of the vector £ and the matrices Q and P in Eqs. < |43[ ) 
and (45 I that determine the Fourier coefficients d„ of t\ (8) 



according to Eq. ( 47 1 



The Fourier decomposition ( 23 1 leads to linear equations 



on the coefficients ao and a n , n> 1, 



R 2 



2K, 



00' 



an I 1 H — 

2dD 2 \ kR 



= <*t|i>-£ 



(0 



kR 



°° rr7/T (e) 

' ' t-i V 77,777. 



fciZ 



<*i|V^> + 



RK, 



00 
■/>() 



(25) 
(26) 



First, one can solve the set of linear equations (26), in- 
dependently of Eq. ( 25 ) , by writing 



where 



M = —K^ 
U n = (h\V n ) + 



R 



-K. 



(e) 



dkD 2 " n0 



(27) 

(28) 
(29) 



for all n,m > 1. Formally, the solution of this system of 
equations is 



(I + M) *[/ 



(30) 



where (I) n ,m — $n,m stands for the identity matrix. Note 
that t\{9) and thus U n are still unknown. 
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Second, substituting the B.C. (15) into the diffusion 



equation ( 12 ) leads to 
D 



C. Exact solution 



lAeh(0) = -l+'-^ . (31) 

R 2 K Or \r=(R,B) 



Expanding the function ip(9) — g t {9) on the complete 
basis of functions {V n (9) — V n (e)} n >i, 



The substitution of Eq. pi into this relation yields ^(°) = dA°) + Yl d n{Vn(9)-V n (e)} (e < 6 < tt), (41) 



-A e t 1 (6) = -u 2 T + u, 2 J2 a n-j^ V n {6), (32) 

n— 1 

where we defined 



CJ = Ry/X/D u 

1 1 

Y + - 
A a 



T 



(33) 
(34) 



a being the inverse of the mean re-adsorption time on 
the surface after a desorption event, as defined in |23| 



dkD 2 
R 



(35) 



The solution of Eq. (32) which satisfies the B.C. 
*i(e)=0,is 

2rp ra s nV n (9)-V n {e) 

ti{6) = u Tg e {9) -u 2^a n — (36) 



kR 



Pn 



for e < 9 < 7r, and t 1 (9) = otherwise (0 < 9 < e). Here, 
R 2 g e (9)/Di is the well-known MFPT for a surface search 
(i.e. when A = 0) [28]: 



g e (9) = In 



1 - cos(0) 
1 — cos(e) 



(37) 



for e < 9 < 7r, and 3 C (6') = otherwise. One can easily 
check that —Agg e = —1 and g t (e) = 0. Using the formal 
expression ( 30 ) for a n and introducing the dimensionless 
function 



h(9) 

w 2 T' 



(38) 



one can represent Eq. ( 36 ) as an integral equation on 



i>(9)=9e(9) + £ 



V n (9)-V n {e) 



Pn 



(39) 



for e < 9 < n, and V'(^) = otherwise. Here, we intro- 
duced 



n 

kR 



(i + apt 1 



(40) 



one obtains a set of linear equations for the unknown 
coefficients {d n } n >i 



n=l 



J2 d n {V n (9)-V n (e)} (42) 

oo / oo \ 



n=l 



1=1 



where we defined the vectors U and £ by their n-th co- 
ordinates (n > 1) 



1 oo / 

P ™ m=l V 
£n = Pn(5£|Kx) e , 



/T (e) P 



p m ' dkD 2 u 2 T 



(43) 
(44) 



and the matrices Q and P e ) by their n-th row and i-th 
column (n, I > 1) 



Qnl — — J]] X nm I m[ , 



m— 1 



J« = (V m {9)\V l {9)-V l {e)) e . 



(45) 



(46) 



As Eq. (42) is satisfied for all 9 € [e, 7r], the coefficients 



d n are found as 

d„ = 



uj 2 (I - lo 2 Q) 1 P 



(47) 



Combining this relation with Eqs. rt37 38 41 ), one finally 



obtains an exact representation for the MFPT 

ti(0) = w 2 P ^ e (0) + - K„(e)}J (48) 

for e < 6» < tt, and t x {6) = for < 9 < e. 

Averaging t\(9) over the whole surface and using the 
relation 

(P„(cos60-P„(cose)|l) e = Pn (P n (cos9)\g e (9)) e = £„, 
we obtain the exact formula for the surface GMFPT: 



(ti) - w 2 T ( 5e ) + £ , 



where (<? e ) is computed by integrating Eq. (37) 

2 \ 1 + cos e 



(Se> - In 



1 — cos e 



(49) 



(50) 
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Equations (48) and (49) are among the main results of 



the paper, and several comments arc in order: 

(i) As expected, in both limits A = and k = oo, we 
retrieve the limit of the 2 MFPTs for the surface search 
process alone: t\(9) — > j^g t {9). 

(ii) A physical interpretation of Eq. (49) is that 



the surface GMFPT is the product of the mean time 
T = A -1 + a -1 for an elementary cycle composed of one 
surface exploration and one bulk excursion, by the mean 
number of cycles before reaching the target. 

(iii) A numerical implementation of the exact solutions 
in Eqs. ( 48 ) and ( 49 ) requires a truncation of the infinite- 
dimensional matrix Q to a finite size N x N . After a 
direct nume rica l inversion of the trunca ted matrices (I 

uj 2 Q) in Eq. 



M) in Eq. fl30|) an d ( 

from Eqs 



(48 



(47), the MFPTs 



49 ) are approximated by truncated series 
(with N termsy _ At a fixed tolerance threshold, higher 
values of A require higher values of N. In spite of the 
truncation, we will refer to the results obtained by this 
numerical procedure as exact solutions, as their accuracy 
can be arbitrarily improved by increasing the truncation 
size N. 



le series are truncated at suffi- 



(iv) The expression (49) is valid for arbitrary target 
size e, provided that tl 
ciently high N. 

Substituting Eqs. 



(38) and (41 ) into Eq. (29), we 



deduce the Fourier coefficients a n of t 2 (r,9) for n > 1 




(51) 



RK 



(<0 



dkDo 



while ao is found from Eq. ( 25 ) . The coefficients a n 



determine an exact representation (23) of the MFPT 
t2{r,9). Therefore one gets a complete exact solution 
of the problem for any starting point. In particular, the 
bulk GMFPT (t 2 ) averaged over uniformly distributed 
starting points in the bulk, reads 



<^2 



2tt 



47ri? 3 /3 



It 7T 

J drr 2 J d6 sin 9 t 2 (r, 



9) =a - 



m 2 

lOdD-2 



(52) 

since the other terms from Eq. (23) vanish due to the 



orthogonality of V n (9). Substituting an expression for 
a>o, one gets 



{'2 




kR 



(53) 

Note that the coefficients a n in Ref. [22] were simply 
proportional to the coefficients d ni up to the third order 
in e. The mixed boundary condition ( 15 ) results in the 



more sophisticated expression (51) for a 



D. Existence of an optimum 

Despite the prediction of [23] that the bulk and surface 
GMFPTs are monotonic functions of A, the exact solu- 
tions prove to admit a minimum, as seen in Fig. [3]d on 
the example of the surface GMFPT. In this section we 
focus on the surface GMFPT and we determine sufficient 
conditions for this GMFPT to be an optimizable function 
of A, which are set by two requirements: (i) desorption 
events should decrease the search time for small enough 
values of A, i.e. (ti) < (ti)x=o = t s ; (ii) the mean surface 
search time is lower than the mean bulk search time, i.e. 

t s < t b . 

The first condition is fulfilled when the derivative of the 
surface GMFPT is negative at A = 0. We first rewrite 
the coefficients U. n as 



U„ 



1 



dkR + A|| D 2 



where 



z n 



Pn n Pm 
rn—1 

oo 

— V X K {€) 



(54) 

(55) 
(56) 



The derivative of the GMFPT at A = is 
'd(t 1 )\ R 4 ( Di 



OX 



a=o D \ \dkRD 2 



(9e) + (£W) )+(t-Z) 

(57) 



This derivative is negative provided that 

J_ (.9.) + (± W) 
Dj low dkR {i-Z) ■ 



D2 > (D2 



( 58 ) 

where we used the inequality (£■ Z) < 0. 

The second bound (ii) is obtained when the surface 
search time R 2 (g e ) /D\ (at zero desorption rate) is lower 
than the search time at infinite desorption rate. Using 
the first-order asymptotic expansion of Eqs. ([2j [3]), this 
condition is explicitly given up to the second order in 
e -C 1 as 



Ei < ( £i 



(59) 



Eh 
Di 



3e(l - elne)(21n(2/e) - 1) 



0(e) 



Combining the above inequalities, one gets a sufficient 
condition for the surface GMFPT to be optimizable for 
e<l: 



Lh 



D 2 I ' D 2 



low 



(60) 



Figure [4] displays the regime of parameters for which the 
surface GMFPT is optimizable. 
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E. Perturbative solution 

The first terms of a perturbative expansion with re- 
spect to e of Eq. (49) can easily be obtained. At first 
order in e, one has 



U n = U<® + 0(e) 



V2n+ 1 



kR 



Q 



» 2 (« + i) 2 i + 4 
i — 

1 kR 



0(e) 



- O (0) + O(e) - : kH S 

ran — ^ ran ' ^ W ~ i M \ , , n u m,n 

n(n + 1)1 + m 



0(e), 



from which 

d n =^[(I-^Q^U^] n + 0(e) 



H{n + 1) «(« + l)+- 2 (lT^) 
One finds therefore 
ip(6) = -21ne + 21n(2sin(6»/2)) 



0(e) 



2 f f m \ 2n + l l-P„(cos( 



Si 



f 1 \ 1 + mJ n(n + l) ?1 ( n + i) +w 2 
+ 0(e). (61) 
Averaging the latter relation over # yields: 



(tl) = w^T 21n(2/e) - 1 



2n + 1 



fei? 



S n ( n+ i) n ( n+ i)(i + 1 ^) +I 



,2 " 



O(e) 



(62) 



Notice that this expression is identical to the pertur- 
bative development at e < 1 for the case of a uniformly 
semi- reflecting sphere, including the target |22) . In the 
limit of very small e, the target is mainly reached from 
the adsorbed (surface) state, while its reactivity from the 
bulk is expected to be negligible. As seen on Fig. [3j the 
smaller the desorption rate A, the larger the domain of 
applicability in e of the perturbative development. Note 
also that in Fig. [3]d, the first-order expression exhibits a 
minimum with respect to the desorption rate. 



IV. COMPARISON WITH THE 
COARSE-GRAINED METHOD 

A. Coarse-grained approach 

In the coarse-grained approach of [53], the bulk, the 
target and the rest of the surface are considered as effec- 
tive states, denoted by b, s and 0, with no inner spatial 
degrees of freedom. An effective set of kinetic equations 
combines four first-order reaction rates: (i) the rate A 



(denoted as /? in [53]) which is associated with the desorp- 
tion (s — > b); (ii) the rate a defined in Eq. (35) as the in- 



verse of the mean re-adsorption time on the sphere, which 
is associated with the effective adsorption (b — > s); (iii) 
the rate k s defined as the inverse of the surface GMFPT 
of Eq. ([2| for the surface search alone (s — > 0): 



.4 



(63) 



where the prefactor A from Eq. ([T]) accounts for the 
difference between the averages over the initial position 
mentioned above; (iv) the rate kb defined as the inverse 
of the bulk GMFPT of Eq. Q for the bulk search alone 
(6 ->• 0). As in [53], we define 



dD 2 e 
irR 2 ' 



(64) 



which is the inverse of the first-order asymptotics in the 
limit e < 1 of Eq. 

In the following, we focus on the MFPT averaged over 
the sphere surface (t\) and compare our exact expression 
Eq. ( 49 ) to the concise approximate expression derived 



in 123 for the surface GMFPT 



(h) * A 



a + A + kt 



ak s + \kb + kbk s 



(65) 



As found in [23], this expression predicts a monotonic 
behavior for the surface GMFPT as a function of the 
desorption rate A. Similar results can be obtained for 
the MFPT averaged over the sphere volume fa)- 



B. Comparison 

Figure [2] shows that the coarse-grained approach and 
the exact solution are in good agreement for the val- 
ues of the adsorption parameter kR = 6.4 • 10 -4 , 6.4 • 
1CP 3 , 6.4 • 1CP 2 , which are used for an analogous plot 
in 23J. However, as soon as kR is large enough, the ex- 
act and approximate curves are significantly different, as 
illustrated on Fig. § (for kR = 1, 10, 100). Finally, 
and most importantly, the exact solution for the surface 
GMFPT (ti) can exhibit a minimum with respect to A 
(see Fig. [3]b) as opposed to the coarse-grained approach 
which always predicts a monotonic behavior. Similarly, 
the bulk GMFPT (t 2 ) from Eq. (pi) also exhibits a min- 



imum with respect to A for the parameters used in Fig. 
[3]d (not shown). 

Note that in all considered cases of Figs. [2] [3] the 
coarse-grained approach underestimates the search time 
at large A. This discrepancy is due to the approximate ex- 
pression of tb used in [53] to estimate kb- As expected, in 
a similar two-dimensional problem, our solution at large 
A coincides with the earlier exact result of [24] for the 
mean exit time when the surface is perfectly reflecting 
(see Appendix [E| . 




FIG. 2: (Color online) Comparison between the coarse- 
grained (symbols) and exact (lines) surface GMFPT (ti) as 
a function of the desorption rate A, for e = 0.02, N = 3 • 10 4 
and several values of kR = 6.4 ■ 10~ 4 , 6.4 • 10~ 3 , 6.4 • 10" 2 
(from [13), with D 2 = Di = 1 (a) and D 2 = 4£>i = 4 (b) in 
arbitrary units in which R = 1. 



FIG. 3: (Color online). Comparison between the coarse- 
grained (symbols) and exact (lines) surface GMFPT (ti) as 
a function of the desorption rate A, for e = 0.02, N = 3 • 10 4 
and several values of kR — 1, 10, 100, with D2 = Di = 1 (a) 
and D2 = 4Di = 4 (b). The dotted curves illustrate the 
perturbative solution ( |62[ ) which is accurate for moderate A 
but strongly deviates for very large A. 



C. Discussion 

We now discuss quantitatively the validity domain of 
the coarse-grained approach and explain qualitatively 
why it fails to reproduce the minimum of the surface 
GMFPTs with respect to A. 

As suggested in [23] , the validity of the coarse-grained 
approach requires the equilibration times for homoge- 
nization within each state to be faster than the other 
time scales of the process. In particular, the equilibra- 
tion rate should be larger than the following: 
(i) the target encounter rate, i.e., 



max(fc s , kb) <C mm (D 1 /R 2 ,D 2 /R 2 ) 1 



(66) 



which implies e< 1. The coarse-grained approach can 
only describe narrow escape situations e <C 1, as men- 
tioned in [35], whereas the exact solution presented here 
is valid for arbitrary e. 

(ii) the inverse of the mean time for re-adsorption on the 
surface after desorption, i.e., 



a = 3kD 2 /R < min(Di/ J R 2 , D 2 /R 2 ), 



(67) 



which implies in particular that kR <C 1. Indeed, as seen 
in Appendix [C] in the regime kR 1 , the spatial cor- 
relations between the starting and ending points of bulk 
excursions are negligible. The condition (67) is satisfied 
in the situations displayed in Fig. [2] 

However, Eq. ([67]) is not longer satisfied in both situ- 



ations of Fig. [3j and one notices that the coarse-grained 
solution does not match with the exact result. This is 
particularly visible in Fig. |3j), where the surface GMFPT 
(ti) exhibits a minimum with respect to the desorption 
rate A. 



The condition (67) for the applicability of the coarse- 



grained approach turns out to be incompatible with the 
existence of the minima of the bulk and surface GMFPTs 
with respect to A. Such a minimum can be attributed to 
the fact that bulk excursions reduce the time loss due to 
the recurrence of surface Brownian motion by bringing 
the particle, through the bulk, to unvisited regions of the 
surface. In the coarse-grained approach, a bulk excursion 
s — > b — > s "consumes" time but brings the particle back 
to its initial effective state s. The assumption that the 
starting and ending points of bulk excursions belong to 
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bulk excursions may be too favorable 




surface diffusion 
is preferred 



(a) 



10"' 



10" 



kR 



10" 



10' 



7p 

6 

5 



Q 4 



bulk excursions may 
be too favorable 



3 
2 





optimizable 
region 



surface diffusion is preferred 



0.05 



0.1 

e 



0.15 



0.2 



FIG. 4: (Color online). The regions of optimality for the 
surface GMFPT (ti). (a) The critical ratio £> 2 /£>i as a func- 
tion of the adsorption coefficient k for a fixed value of the 
target half-width e = 0.02; (b) The critical ratio D2/D1 as 
a function of the target size e, for the adsorption coefficient 
kR — 10. Below the lower bound (dashed red line), surface 
diffusion is preferred. Above the upper bound (solid blue 
line), the GMFPT is smaller in the desorbed state than in 
the adsorbed state. In between, the surface GMFPT is an 
optimizable function of A. Series are truncated at N = 10 4 . 



tion is presented in Appendix [B] in which the agreement 
is found within the domain of applicability of the pertur- 



bative expansion (introduced in Sec. HIE). 



V. CONCLUSION 

We have obtained an exact expression for the MFPT 
and their spatial averages, called bulk and surface GMF- 
PTs, for the process studied in [23]. Compared to [TSUI], 
the introduction of the surface binding rate k allowed one 
to avoid using the ejection distance a (i.e. to set a = 0) 
after each desorption events. In contrast to the statement 
of [23], we have shown that the bulk and surface GMF- 
PTs can be optimized even in this situation and that this 
optimality property is not related to the non-locality of 
the intermittent process considered in [T2 [T3J . 

These exact results can be extended in several direc- 
tions to include the following: (i) the search for a semi- 
reflecting target, with an adsorption parameter different 
from the rest of the surface (see Appendix [D| ; (ii) the 
2D search for an angular aperture on the boundary of 
a disk (see Appendix |E| ; and (iii) the biased search for 
an arbitrary number of regularly spaced targets over an 
otherwise semi- reflecting annulus (2D) or cone (3D), fol- 
lowing the method of [22] . Notably, even for a 3D search 
for a purely reflecting target and for a 2D search with 
a bulk diffusion coefficient D2 smaller than the surface 
diffusion coefficient D\, the surface GMFPT can be an 
optimizable function with respect to the desorption rate 
(see Figs. [6] and [7] below). 
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Appendix A: Boundary condition on the MFPT 

The diffusion equation and the boundary condition for 
the occupation probability distribution in the bulk ap- 
pear at first sight to be different in [33] and [22]. In 
this section, we show that these two sets of equations are 
in agreement and lead to the same radiative boundary 
condition (Jl5|) for the MFPT. 

In 



(671 is a c 



lallenging question. A first step in this direc- 



we defined a different set of equations for the 
conditional probability distribution to include a radial 
ejection distance a after each desorption event, in the 
presence of a velocity field. With the shorthand notations 
of Sec. Ill A| the forward advection-diffusion equation of 
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reads: 
dp(r, 



dr 
dp(r, 6 
ds 

+ X 



= -kp(R,9) + 



v(R) 



D 2 
R 



D 2 
v(r) 
D2~r 



P(R,0) 
p(r,6) 



R-a 



S 3 (v-(R-a,9))p(9), 



(Al) 



(A2) 



where v(r) is a radial velocity field positive for an out- 
ward drift. The other diffusion equations on the condi- 
tional probabilities are the same as in |22j . 

Notice that Eq. (|A1| does not involve the desorption 



rate A. However there is no contradiction with Eq. (16 1, 



as Eq. (Al) (without drift) leads the same backward 
boundary Eq. ( 19 ) as long the appropriate limit for a = 



in the Dirac function in Eq. ( A2 ) is 

R , 



R 



R-a 



S(r -(R-a, 6))p(r, 9)r 2 dr ^—^ R 2 p(R, ( 



It can be proved that this condition is required from nor- 
malization of the probability density. 



Appendix B: Analytical agreement between the 
coarse-grained and exact solutions 

Figure [2] shows a good agreement between the exact 
and coarse-grained expressions (49 65 ) . However, finding 



an explicit analytical relation between these two expres- 



sions under the general conditions (66 67) seems non- 
trivial. 

We focus here on the following specific successive lim- 
its: (i) small target extension (e <C 1); (ii) low-desorption 
rate regime (A <C k s ); (hi) and intermediate range for the 
adsorption rate: kb <C a <C min(_D 2 /-R 2 , Di/R 2 )- 

On one hand, the coarse-grained expression in these 
limits reads 



t s + 0[e 



(Bl) 



where we have used that kb/k s <C 1 and kb/a ^ 1 for a 
fixed value of a at sufficiently small t< 1. 

On the other hand, in the limit e <C 1 the perturba- 
tive expansion of Eq. (62) is accurate. In the above- 



mcntionncd limits it reads 



A 



1+- )t s + 0(e,X/k s ,kR), 



(B2) 



where we have used that 



E 



2n + 1 



k-R 



,tt n(n + l) n(n + l)(i + ^) 
« l + 0(kR,R 2 X/D 1 ) 1 



+ 



\R 2 n 
Dx kR 



(B3) 



and that the condition R 2 X/Di <C 1 is guaranteed from 
the condition X/k s <C 1. 

The identification of the first order terms in Eqs 



(Bl 



B2 ) shows analytically the agreement between the coarse- 



grained and exact solutions within a range of parameters 
which necessarily satisfies Eqs. (66 67). Notice that the 



argument presented here relies on the condition k^ /a <C 1 
which is not satisfied in the situations represented on Fig. 
® 



Appendix C: Measure of correlations 

In this section we quantify the spatial correlations be- 
tween the starting and ending points of a bulk excursion. 
We then provide the range of values for k in which the 
spatial correlations are negligible. The probability den- 
sity n(6*) for a particle initially started from the surface 
state at the angle 9o — 0, </>o = to first return on the 
surface to any point (6,4>) (with </> € [0, 2ir}) is [22] 



n(0) 



sin 9 



1 



2n+ 1 



P„(cos( 



(CI) 



kR 



The cumulative probability distribution for the reloca- 
tion angle 9 E [0, 7r] is the integral of the probability 
density 



F(0) = 



1 



n(6»')sin6i / d9'. 



(C2) 



By analogy with the Kolmogorov-Smirnov test |25j . we 
propose to measure the spatial correlations between the 
starting and ending points of a bulk excursion by the 
norm 



N K = max \5F(9)\ 

06[O,tt] 



(C3) 



with 



5F(e) = F(e)-F a (9) = -Y, 



1 v 2 ^ P„_i(cos6>) - P n+ i(cos( 



n=l 



k:R 



(C4) 



and F u (9) = (1 — cos6')/2 is the cumulative distribution 
for uncorrelated random relocation on the sphere. This 
leads a correlation angle O which is defined as the solu- 
tion of the equation SF(Q) — Nk- 

As shown on Fig. [5] the spatial correlation is negligible 
(Nk -C 1,6 « 1) as long as kR < 1. In particular for 
the reference values kR = 6.4 • 10~ 4 , 6.4 • 10~ 3 , 6.4 • 10~ 2 
used in [23], the norm Nx is smaller than 0.05 and the 
correlation length is nearly constant at O w 1.2. 



Appendix D: Generalization to semi-reflecting 
targets 

In this section, we briefly generalize our method to the 
case of a semi-reflecting target, for which the B.C. reads 
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FIG. 5: (Color online). The norm Nk (blue solid line and 
circles) and the correlation angle O in radians (red dashed line 
and squares) as functions of the adsorption coefficient k, for 
the truncation size N = 3000 (lines) and N = 1000 (symbols). 



t 2 (R,0) = 



k t dr \r=(R,e) 

k dr |r=(H,i 



(0 < 6 < c), 
(e < 6> < tt). 



(Dl) 



This general description includes the following cases: 



• kt — oo is the case considered so far in the present 



article of a fully adsorbing target, Eq. ( 15 ) 
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FIG. 6: (Color online). The surface GMFPT (ti) as a func- 
tion of the desorption rate A, for several values of the target 
adsorption parameter kt (for particles hitting the target from 
the bulk), with e = 0.02 and D2 = Di = 1 in units in which 
R — 1. The non-reactive part of the surface is semi-reflecting 
with an adsorption parameter k = 100. Note that the tar- 
get is perfectly reactive for particles adsorbed on the surface. 
Series are truncated at iV = 3 • 10 4 . 



• kt — k is considered in Ref. [22] : 



• kt = corresponds to a target which is fully reflect- 
ing for particles hitting the target from the bulk. 
Using the mixed boundary condition (Dl I, the projec- 
tion of a series representation ( 23 1 for t^R, 0) onto the 
basis {V n (6)} n > becomes 



ti(0)V n (0)d/i(0)- 



1 dh 

dr \r=(R,6 



V n (8)dfi(6) 



(D2) 



1 \ dU 



1 



k k t J dr \r 



V n {9)dfi(9), 



which replaces Eq. (24), with dfi{6) = ^sin8d9. This 



leads to the following equations on coefficients a n : 




From these equations, we extend the definition of M„ 
and U m from Eq. (28 1 to 



Following the same steps as in Sec. IIIB| one gets an 
integral equation on the dimensionless MFPT ip(9): 



M nm = 6^^ + m (±-±)K&, (D5) 



K ( :l (D6) 



V n (6)-V n (e) 



Xr, 



Pr> 



(D7) 



( R 2 



\dD 2 T \kR 



12 



which generalizes Eq. ( 39 1 . Expanding this function onto 
the basis {V n (6) 



1 °° ( 

rn—1 \ 



V n (ejJ yields Eq. {47}, with 
R 2 



1 



1 



K 



(<0' 



(D8) 

which generalizes Eq. (43). This relation can also be 
written as 



U n — Z„ 



Di 



dD 2 (l + X/a) 



1 

kR 



1 



(D9) 



where Z n and W n are still defined through Eqs. (55l 56 1. 
Other quantities and representations remain unchanged. 
Repeating the computation of the derivative of (ti) at 
A = 0, one gets the lower bound as 



D2 
D 1 



1 (ge) + (S-W)(l-k/k t ) 



dkR 



(D10) 



which extends Eq. (58). 

One can see that the change in the boundary condi- 
tion, i.e., extension from Eq. (15) to Eq. (Dl), does not 



affect the method and the structure of the solution. As a 
consequence, the conclusions on the optimility of the sur- 
face GMFPT remain qualitatively unchanged, although 
values of the lower bound may be different. Note that the 
determination of the upper bound requires the expression 
of the surface GMFPT for a semi- reflecting target in an 
otherwise reflecting sphere, which is still unknown. 

As shown on Fig. |6j optimization in A remains possible 
even in the case of a target which is fully reflecting for 
particles hitting the target from the bulk. 



Appendix E: Application to a disk 



for the Laplace operators: 
A 



dr 2 



ld_ 
r dr ' 



A/, 



dr, 



The eigenfunctions V n (9) and eigenvalues p n of the an- 
gular Laplace operator Ag become 



V n {9) = V2 cos n9, p n 



(n > 0). (El) 



The inner scalar product (21 ) is replaced by 

(f,g)^{f\g)e = - r f(0)g(0)d9. (E2) 



In Table [TTJ we summarize the expressions for the matrix 
elements K~mn and Imn which replace those from Table 



[IJ One has also to replace Eq. ( 37 1 by 

gM = l(0-e)(2n-e-9). 



(E3) 



1 

0.9 

^ 0.7 

0.6 
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1.09-10" 2 




-■-k = 


1.09-10" 1 




~t~ k = 


1.09 



10 



10" 



10' 

A7k 



10' 



10° 



10' 



V n {6) 


V2cos(n9) 






<5e|l}e = <5e>e 


hi---? 


£» (n > 1) 


— — { (i" — e) cos(ne) + sin(ne) / n} 


-"■oo 


(tt — 


K% (n>l) 


-^sin(ne) 


Kit (n > 1) 


1 / sin(2ne) \ 

- 7T — € 7i 

7r \ 2n / 


K^m n,m> / n 


1 ( sin((m — n)e) , sin((m+n)e) \ 
7r \ m — ?i 1 m + n / 


Ini (n > 1) 


1 ( ^ r i sin 2ne \ 

- 7T — e H = — ) 


Inm 771,71 > 1,771 7^ 71 


^sinfm^ sinW) 

2 cos(ne) --cos(me) 2 

— — 75 t — 771 

7T n^—m,^ 



TABLE II: Summary for the 2D case of the quantities in- 
volved in the computation of the vector £ and the matrices Q 



and M in Eqs. (43 45 I that determine the Fourier coefficients 



d n of t\(9) according to Eq. ( |47| ). 

In 2D, the diffusion equations on the MFPT are iden- 




FIG. 7: (Color online). The surface GMFPT (ti) in the disk 
as a function of the desorption rate A, for e = 0.02, D2/-D1 = 
1, and (a) kR = 1.09- 10" 2 , 1.09- 10"\ 1.09 (corresponding 
to a = O.lfc;,, kf,, lOfc;,) and (b) kR = 1, 10, 100. Series are 
truncated at TV = 3 ■ 10 4 . 



The method of derivation and the remaining quanti- 
ties are not modified. In particular, Eq. (48) for the 
tical to Eqs. (Tl2j fl3|) provided the change of definition MFPT ti(9), Eq. Q for the surface GMFPT (tt) and 
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bulk excursions may be too favorable 



1.5 



surface diffusion 
0-5 is preferred 



optimizable region 



Eq. ( 58 ) for the lower bound on the diffusion coefficient 
ratio (-D2/-Di)i are applicable for the disk. In turn, the 
asymptotic relations on the exit time (pi [31 are modified: 



t s — (ti)x=o — 



R 2 (n-e) 3 



D x 3?r 
j R?_ 
" D 2 



ln(2/e) + Oil) 



(E4) 
(E5) 



10" 4 10" 3 10" 2 10" 1 10° 10 1 10 2 10 3 

k 

FIG. 8: (Color online). The regions of optimality for the sur- 
face GMFPT (ti) for a disk of radius R = 1 with an aperture 
of half-width e = 0.02. Below the lower bound (solid red line) 
surface diffusion is preferred. Above the upper bound (dashed 
green line), the surface GMFPT is higher in the adsorbed 
state than in the desorbed state. In between, the surface 
GMFPT is an optimizable function of A. Series are truncated 
at N = 10 4 . 



that leads to the following expression for the upper bound 
on the diffusion coefficients 

Figure [7] illustrates the fact that the surface GMFPT 
is an optimizable function of the desorption rate A in the 
range of parameters represented on Fig. [8] 
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